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Abstract. 

We propose a jet model for the low/hard state of Galactic black-hole X-ray sources that can explain the energy 
spectra from radio to X-rays. The model assumes that i) there is a magnetic field along the axis of the jet; ii) the 
electron density in the jet drops inversely proportional to distance; and iii) the electrons in the jet follow a power 
law distribution function. We have performed Monte Carlo simulations of Compton upscattering of soft photons 
from the accretion disk and have found power-law high-energy spectra with photon-number index in the range 
1.5 - 2 and cutoff at a few hundred keV. The spectrum at long wavelengths (radio, infrared, optical) is modeled 
to come from synchrotron radiation of the energetic electrons in the jet. We find flat to inverted radio spectra 
that extend from the radio up to about the optical band. For magnetic field strengths of the order of 10^ — 10^ G 
at the base of the jet, the calculated spectra agree well in slope and flux with the observations. Our model has 
the advantage over other existing models that it also explains many of the existing timing data such as the time 
lag spectra, the hardening of the power density spectra and the narrowing of the autocorrelation function with 
increasing photon energy. 
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K>" 1. Introduction 



, Galactic black-hole binaries are classified by their X-ray 
■ features such as the strength and temperature of the soft 
multi color black-body component (Mitsuda et al. 1984), 
the hard X-ray emission, the X-ray luminosity and the 
timing properties. Several states have been identified to 
characterize black-hole accretion (for reviews see Tanaka 
& Lewin 1995; van der Klis 1995; Nowak 1995; Poutanen 
1998; McClintock & Remillard 2005). A source is in the 
so-called low/hard state when it exhibits low X-ray lumi- 
nosity (less than a few % the Eddington luminosity Le), 
weak or absent thermal component, dominant hard X-ray 
emission with a cutoff at ~ a few hundred keV and high 
rms variability. 

The soft thermal component is widely accepted to 
come from an optically thick, geometrically thin accretion 
disk (Shakura & Sunyaev 1973), while the hard X-ray tail 
is usually assumed to come from inverse Compton scat- 
tering of soft photons off energetic electrons (e.g. Sunyaev 
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& Titarchuk 1980). The location of the energetic electrons 
and, more generally, the accretion geometry are poorly un- 
derstood. One possibility is that the disk extends all the 
way down to the last stable circular orbit and that there is 
a hot rarefied "corona" that lies above and below the thin 
disk (Galeev et al. 1979; Haardt & Maraschi 1993; Haardt 
et al. 1994; Stern et al. 1995; Poutanen & Fabian 1999). 
It is also possible that the disk is truncated at a larger 
radius and the inner region is filled with optically thin, 
geometrically thick, two-temperature plasma (Shapiro et 
al. 1976; Ichimaru 1977; Rees et al. 1982; Narayan & Yi 
1994; Abramowicz et al. 1995; Esin et al. 1998). 

During recent years, evidence has mounted that, when 
in the low/hard state, X-ray binaries exhibit steady core 
radio emission (Hjellming & Han 1995; Mirabel et al. 1998; 
Fender 2001; Gallo et al. 2003). Moreover, in some cases 
a jet-like structure has been resolved (Mirabel et al. 1992; 
Dhawan et al. 2000; Stirling et al. 2001; Martf et al. 2002; 
Fuchs et al. 2003). The spectrum in the radio is flat to 
inverted and seems to extend to the infrared or even be- 
yond (see, e.g., Hannikainen et al. 1998; Fender et al. 2000; 
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Corbel et al. 2000). Because of the high brightness tem- 
peratures, non-thermal spectra and, in some cases, high 
degree of polarization, radio emission is believed to come 
from synchrotron radiation of relativistic electrons in a jet. 
Furthermore, an interesting radio/X-ray correlation has 
been established to hold through simultaneous radio and 
X-ray observations (Gallo et al. 2003) which can be ex- 
tended to include the mass of the black hole and to define 
the "fundamental plane" of black hole activity (Merloni 
et al. 2003; sec also Falcke et al. 2004). 

Inverse Compton scattering by relativistic electrons in 
a jet has been proposed as a mechanism for the production 
of X-rays and 7-rays in X-ray binaries (Band & Grindlay 
1986; Levinson & Blandford 1996; Georganopoulos et al. 
2002; Romero et al. 2002). The possibility that optically 
thin synchrotron emission from the jet results in the hard 
X-ray tail has also been explored (Markoff et al. 2001; 
Vadawale et al. 2001; Corbel & Fender 2002; Markoff et 
al. 2003; sec also Markoff & Nowak 2004 for constraints 
on jet models via reflection). 

Recently, Reig et al. (2003) (hereafter Paper I) pro- 
posed a jet model that can explain the X-ray energy spec- 
tra and the dependence of time lags on Fourier frequency 
(see for example Miyamoto et al. 1988; Nowak et al. 1999; 
Ford et al. 1999) in terms of inverse Compton scatter- 
ing of soft, disk photons by energetic electrons in the jet. 
The density of the electrons in this simple jet model is 
assumed to drop inversely proportional to the vertical dis- 
tance z from the black hole. Giannios et al. (2004) (here- 
after Paper II; see also Kylafis et al. 2004) further assumed 
that the electrons close to the core of the jet are more en- 
ergetic than those at its periphery and showed that both 
the hardening of the high-frequency power spectra and 
the narrowing of the auto-correlation function with pho- 
ton energy observed in Cygnus X-1 (Nowak et al. 1999; 
Revnivtsev et al. 2000; Maccarone et al. 2000) are repro- 
duced. 

Although successfiil in explaining a large number of 
spectral and timing properties in the X-ray domain. 
Papers I, II did not deal with the part of the electromag- 
netic spectrum that revealed the existence of the jet in the 
first place, i.e. the radio. Radio observations on the other 
hand contain valuable information (spectral slope, flux) 
that can place constrains on any jet model. In this work, 
we model the radio emission in terms of synchrotron ra- 
diation of the electrons in the jet. The main modification 
with respect to the jet model used in Papers I, II is that we 
take the electrons to have a power law energy distribution 
function. We first make sure that this modification does 
not alter any of the previously derived results in the X-ray 
domain. This puts significant constraints on the index of 
the power law energy distribution of the electrons. Then 
we show that flat to inverted radio spectra are naturally 
reproduced by the model. Finally, we apply the model to 
the broad band spectra of XTE J1118-I-480 and Cygnus 
X-1. 



2. The model 

2.1. Description of the Jet 

The jet model that we will use here is built on the model 
of Paper I. Here we describe the characteristics of the jet. 
We assume that the jet is accelerated close to its launching 
region and that it has constant velocity iiy (say along the 
2-axis). Furthermore, we assume that the electron density 
in the jet drops inversely proportional with distance z 
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where no is the density and zq is the height at the base of 
the jet respectively. Mass conservation M oc f |[n(z)_R^(z) 
then determines the polar radius R of the outer edge of 
the jet as a function of z 
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where i?o is the radius at the base of the jet. The Thomson 
optical depth along the axis of the jet is given by the 
integration of dr = n{z)aTdz from zo to H (where H is 
the extent of the jet) 



Til = noCTT ^0 In (— )■ 



(3) 



In the previous works (Papers I, II) the magnetic fleld 
had been assumed, for simplicity, to be parallel to the axis 
of the jet (z-axis) and homogeneous. Here, we will keep 
B II z but take its strength to vary with z as dictated by 
magnetic flux conservation along the jet: B{z)itR^{z) = 
const. So we have for the 2;-dependence of the magnetic 
field 



B{z)=Bo^. 

z 



(4) 



The main modification with respect to Papers I and 
II is in the energy distribution of the electrons. In Paper 
I, the electrons had been taken monoenergetic with their 
velocities to have a constant spiraling perpendicular com- 
ponent v±. In Paper II, v± had been assumed to drop 
linearly with polar distance by a factor of a few (resulting 
in a jet with a core "hotter" than its periphery), explain- 
ing the hardening of the high-frequency power spectrum 
with increasing photon energy in Cygnus X-1 (Nowak et 
al. 1999). 

Here we will assume that the electrons have a distri- 
bution of v± that extends to ultra-relativistic electron ve- 
locities. We will describe the distribution in terms of the 

Lorentz factor 



j = l/^l-{vl+vl)/c-^ 



(5) 



and will assume -as it is customary in many jet models 
(e.g. Blandford & Konigl 1979)- a power law form 



N{j)d-f oc 7 °'d'-f. 



(6) 



The distribution extends from 7inin to 7max with 7max 
7min. In total, three parameters are needed to determine 
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the electron distribution inside the jet i.e. the index a, 
7inin, and 7max- Since the distribution itself is expressed 
in terms of the Lorentz factor 7, we prefer to keep 7min as 
a parameter instead of v± (of course these quantities are 
related through Eq. |SJ|). 

With this model at hand, the task is to explore its 
spectral properties in essentially the whole electromag- 
netic spectrum so as to compare with those observed in 
black-hole candidates in the low/hard state. Synchrotron 
emission from the energetic electrons inside the magnetic 
field given by Eq. @ will be shown to dominate from the 
radio up to about the optical wavelengths, while we model 
the hard X-ray part of the spectrum as a result of inverse 
Compton scattering of a soft-photon input at the base of 
the jet by the same population of electrons. We have used 
the Monte Carlo method (Pozdnyakov et al. 1983) to sim- 
ulate inverse Compton scattering. 

2.2. Simulations of inverse Compton scattering in the 
Jet 

Our code is very similar to the one described in Papers I 
and II. Here, we describe only the emission of the soft pho- 
tons. Multi-color black-body photons (i.e. accretion disk 
photons) of maximum temperature Tbb are assumed to be 
emitted from the underlying disk with an upward isotropic 
distribution and temperature that scales with polar radius 
as T oc (Shakura & Sunyaev 1973). Furthermore, 

synchrotron flux from the jet also contributes to the soft 
photon input to be Comptonized and is taken into account 
when the model is applied to observations (see Sect. 6). 

3. Model parameters 

The present model has several free parameters, most of 
which refer to physical quantities of the jet. These are: 
the temperature Tf,t,, the extent H of the jet, the radius 
Rq and the height zq of the base of the jet, the Thomson 
optical depth r|| along the axis of the jet, the velocity of the 
jet and the strength of the magnetic field at the base of 
the jet Bq. The electron distribution is specified through 
the minimum Lorentz factor 7niin, the maximum Lorentz 
factor 7niax and the exponent of the electron distribution 
a (see Eqs. (5) and (6)). 

The values of the parameters that reproduce quite well 
several properties of black-hole candidates in the low/hard 
state (with special emphasis to Cygnus X-1) are called 
reference values in this work and they are: fcTbb — 0.2 
keV, H = 5 ■ lO^rg, Rq = lOOrg, Zq = bvg, t„ = 2.5, 
W|| = 0.8c, Bo = 3- 10^ G, 7„iin = 2.1, 7^ax - 500 and 
a = 4, where Vg = GM/c^ ~ 1.5 x 10^ cm corresponds to 
the gravitational radius of a 10 solar-mass black hole. 

Since the jet is mildly relativistic, the results also de- 
pend on the angle 6 of observation with respect to the 
jet axis. As in Paper II, we will focus in an intermediate 
range of observing angles 0.2 < cos6' < 0.6. Practically, 
for the Monte Carlo simulation this means that we count 
only photons that leave the jet in this range of angles. 



Despite the fact that the parameter space is rather 
large, all the parameters can be constrained since the 
model is tested against (and succeeds in reproducing) a 
large number of timing and spectral properties of black- 
hole X-ray binaries. It is encouraging that none of the 
parameters needs to be fine tuned in order to match the 
observations. The last statement will be quantified in the 
next sections, where the results are presented. 

3.1. The ejection rates 

Using the reference values, one can calculate the number 
density ng of the electrons (plus possibly positrons) at the 
base of the jet (Eq. (3)). Assuming a neutral jet and / 
electrons per proton (which also means / — 1 positrons), 
the ejection rate in the jet is given by the expression 

M = no^i?2«||(me-H^^^). (7) 

In the case of absence of pairs (i.e. / = 1) and us- 
ing the reference values of the parameters one finds an 
ejection rate M = 8.2 x 10^^ g • sec~^. This value is su- 
per Eddington but not necessarily unrealistic. In the case, 
however, of a pair dominated jet (i.e. / 3> 1) the ejection 
rate can be lower by more than three orders of magnitude 
and much more "reasonable" . At this point, our knowledge 
on the composition of the jet is poor and both baryon and 
pair dominated jets are allowed by observations (see for 
example Fender 2003). 

3.2. The spectral index of the electron distribution 

In this work we treat the power-law index of the electron 
distribution a as a free parameter and we do not address 
the issue of particle acceleration in the jet. However, the 
values of a > 3 (see next Section) that we use in this work 
are unusually large (with respect to what is expected from 
shock acceleration for example). Here, we show that a is 
not the spectral index of the injected electrons but that of 
the (steady-state) emitting electrons after taking into ac- 
count radiative (synchrotron and inverse Compton) cool- 
ing. 

Defining the "cooling Lorentz factor" 7c of the elec- 
trons as the Lorentz factor at which the radiative cooling 
time scale t-cad equals the dynamical time scale ^dyn, where 

, ^ , ^syn^Compton /q\ 

t-dyn — ~ — ^rad — "7 ~J~~7 i \P} 

C t-syn I t'Compton 

one can estimate 7c. The synchrotron cooling time tgyn 
can be calculated straightforwardly with the use of Eq. 
(4), while the Compton cooling time can be calculated 
numerically with the Monte Carlo method (see Sect. 2.2). 
For the reference values of the parameters of the jet it is 
estimated to be 7c ~ a few ~ 7min- So, it is approximately 
correct to assume that the whole electron distribution is 
in the fast cooling regime and can be modeled as a power- 
law. 
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E(keV) 

Fig. 1. Emergent photon- number spectrum multiplied by 
E'^ from our jet model for different values of the minimum 
Lorentz factor 7min of the electron distribution. The solid 
curve corresponds to the reference values of the param- 
eters. The dotted, dash-dotted, dashed and dash-double- 
dotted curves correspond to 7min = 1-8, 2, 2.3, and 2.5 
respectively. The rest of the parameters are kept at their 
reference values. The solid gray curve corresponds to the 
illustrative example where 7min varies with z according to 
Eq. ®. 




E(keV) 

Fig. 2. Emergent photon- number spectrum multiplied by 
from our jet model for different values of the power-law 
index a. The dotted, solid, and dashed curves correspond 
to a = 5, 4, 3 respectively. We also plot (dash-dotted 
curve) the case where a — 2 and 7min — 1.67. The rest of 
the parameters are kept to their reference values. 

4. Results in the X-ray domain 

In this section, we show our results concerning the spectral 
properties of the electromagnetic radiation in the X-rays. 
The temporal properties are discussed very briefly since a 
detailed description has already been given in Paper II. 



4.1. Energy spectra 

In Fig. 1 (solid curve) we plot the energy spectrum of 
the emerging radiation as calculated by our Monte Carlo 
simulation for the reference values of the parameters (see 
Sect. 3) and for a black-body soft photon input ^. Two 
components are evident in this plot. Below ~ 1 keV the 
dominant component is the soft photon input, while for 
E >1 keV the Comptonized component dominates. The 
latter component can be well fitted with a power law of 
photon- number index F ~ 1.7 with an exponential cutoff 
at Ec^t ^ 300 keV. 

In Papers I, II we have explored how the spectrum 
depends on some of the parameters such as the optical 
depth along the jet t\\ or the extent of the base of the jet 
Rq. The result was that increasing r|| or Rq, while keeping 
the rest of the parameters fixed, the spectrum becomes 
harder but with essentially the same energy cutoff. This 
is expected since in both cases there is an increase of the 
mean number of scatterings that the electrons experience 
during their random walk inside the jet and gain more 
energy on average. The cutoff, on the other hand, is mainly 
determined by the energetics of the electrons which does 
not depend on r|| or Rq. The same conclusions have been 
checked to hold in this version of the model. 

Of more interest is to explore how the energetics of the 
electrons influence the emerging spectrum. For a power- 
law distribution of electrons, both 7inin and the index a 
are relevant quantities for this exploration. For a large 
enough, as is the case here (see below), the exact value 
of 7max is not important. It is qualitatively expected that 
by increasing 7,„in or making the distribution flatter (i.e. 
decreasing a), the spectrum becomes harder and the ex- 
ponential cutoff appears at higher energies. 

In Fig. n we plot the emerging spectrum for different 
values of 7min, while the rest of the parameters are kept at 
their reference values. One can clearly see the hardening of 
the spectrum with increasing 7i„in- More quantitatively, by 
increasing 7min in the range 1.8 — 2.5, the photon number 
index F of the hard X-ray slope hardens from ~ 1.9 to 
~ 1.5 which is well within the observed range for black- 
hole binaries in the low/hard state. For simplicity, we have 
assumed that 7inin is constant along the jet. This, however, 
may not be the case and 7i„in may vary as a function of z. 
We have checked that if 7min changes rather slowly with z, 
the resulting spectra do not change much. As an example, 
in Fig. ^ we plot the X-ray spectrum for 

7min = 7i(l + 0.8v/W^), (9) 

where 71 = 1/ y^l — f | /c^. 

^ Since in this Section we explore how the model parame- 
ters affect the X-ray spectra (and mainly the X-ray slope and 
high energy cutoff), the energy distribution of the soft photon 
input is not very important. In Sect. 6, however, where the 
model is applied to data, all relevant soft photon sources such 
as multi-temperature disk and synchrotron photons are taken 
into account. 
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We also show (in Fig. O how the emergent spectrum 
depends on the index a of the electron distribution. As 
a decreases, the spectrum gets harder with a cutoff at 
higher energies. Notice that for a = 2 we use the min- 
imum value one could possibly have for 7niinj namely 
7min = — t;|/c2 = 1.67 (see Eq. (5)). One important 

conclusion that can be drawn from this Figure is that if 
the electron distribution is rather flat (say a ~ 2), the high 
energy cutoff becomes too high (i?cut — 1 MeV or more) in 
comparison with the one typically observed in black-hole 
binaries in the low/hard state (~ a few hundred keV) for 
any choice of 7min- This consideration requires a>3 if this 
model is to agree with spectral observations. As we have 
already discussed in Sect. 3.2., such a steep spectral in- 
dex is the result of fast radiative (synchrotron and inverse 
Compton) cooling of the electrons. 

4.2. Temporal domain 

In Papers I and II, it has been shown that the jet model 
can account for a number of temporal properties of X-ray 
black-hole binaries. Since in this work we have modified 
the energetics of the electrons, it is important to check 
whether the previously derived results are altered or not. 
For example the Fourier frequency dependence of the time 
lags between two different energy bands is found to remain 
essentially the same as that shown in Papers I and II, 
in accordance with observations (see for example Nowak 
et a. 1999; Ford et al. 1999). Furthermore, the width of 
the autocorrelation function has been checked to decrease 
with photon energy (as pointed out observationally by 
Maccarone et al. 2000). Finally, it was shown in Paper 
II that, assuming that more energetic electrons lie in the 
core of the jet in comparison to its periphery, a harden- 
ing of the high-frequency power spectrum with increasing 
photon energy is expected. If we assume -in the context of 
the current model- that the minimum Lorentz factor 7niin 
decreases with polar distance, a similar hardening of the 
power spectrum is derived. 



5. Radio emission from the jet 

The existence of a steady jet when a black-hole binary 
is in the low/hard state has been strongly suggested 
through radio observations (Hjelling & Han 1995; Mirabel 
& Rodriguez 1999; Fender 2001). Furthermore, in some oc- 
casions a jet like structure has even been resolved (Mirabel 
et al. 1992; Dhawan et al. 2000; Stirhng et al. 2001; Marti 
et al. 2002; Fuchs et al. 2003). Our simple jet model has 
succeeded in explaining a large number of spectral and 
timing properties of these sources in the X-ray region in 
terms of inverse Compton scattering of soft photons by 
energetic electrons in the jet. Here we will focus on longer 
wavelengths and explore the predictions of our model in 
this region of the electromagnetic spectrum. 



5 

5.1. The spectral slope in the radio 

The coexistence in the jet of energetic electrons and mag- 
netic fields indicates synchrotron emission as the dominant 
radiative mechanism in the radio. The study of radio emis- 
sion is more straightforward in a frame that is comoving 
with the emitting medium (i.e. the jet). Suppose that an 
electron has a Lorentz factor 7 in the lab frame (i.e. the 
rest frame of the black hole) . Then a Lorentz transforma- 
tion gives the Lorentz factor in the comoving frame 

7co = ^l-vpc^j. (10) 

Since we have assumed that v\\ is constant in the jet, 
7co c< 7. Stated in other words, we still have a power 
law distribution for the electrons in the comoving frame, 
namely 

iV(7co)d7co = C7eo"rf7co. (11) 

For the rest of this work, and for simplicity in the no- 
tation, we will drop the subscript "co" from the Lorentz 
factor. The constant C can be found if we integrate the 
last expression from 7min to 7max ^ 7min and equate it to 
the comoving electron density. Doing so we have 

C^n^l-vyc-^{a- 1)^:^7^, (12) 

where we have also assumed that a > 1, which certainly 
holds, and n is given by Eq. (1). 

The properties of the jet (i.e. density, magnetic field 
strength, energetics of the electrons) may vary only along 
the axis of the jet but not along the perpendicular direc- 
tion. As is customary, we can exploit this symmetry by 
dividing the jet into "slices" across the z-axis and calcu- 
lating the synchrotron emission and absorption at each 
slice separately. Adding the contribution of all slices, we 
get the total emitted power. 

For a relativistic power law distribution of electrons, 
the emitted synchrotron power per unit frequency dP{v) 
in a shce in the range (z, z + dz) is (see, e.g., Eq. (6.36) 
of Rybicki & Lightman 1979) 

dP{v) = C^B'^^'^^n^l - vpc^iy^R^dz, (13) 

where Ca depends only on the index a (defined in Eq. 
(6)), while B = B{z), n = n{z), and R = R{z) are given 
by Eqs. CJ, and ©, respectively. 

The synchrotron radiation is strongly self-absorbed be- 
low a characteristic frequency (the turn-over frequency) 
vt. Synchrotron absorption coefficient is given by the ex- 
pression (see Eq. (6.53) in Rybicki & Lightman 1979) 

a, = A^n-l'^^^B^v-^, (14) 

where depends only on the index a. The turn-over 
frequency in a slice of the jet can be estimated as the 
frequency for which the optical depth to synchrotron ab- 
sorption, across a radius of this slice becomes unity, i.e. 
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a^^R{z) ~ 1. Solving this expression for vt and using Eqs. 
(El, 0), and lO we have 



, + 4 „+4 Ro 

'mm U 



■ + 3 
'+4 



(15) 



In estimating the turn-over frequency by setting 
a^^R{z) ~ 1, we actually assume that the jet thickness 
R{z) is much less than its height z. This does not hold 
close to the base of the jet where z < R{z). From Eq. (2) 
we can verify that z ^ R{z) when z ^ lO'^rg. Since the 
jet extends up iJ = 5 • lO^rg, our assumption holds along 
most of the jet. We will return to this discussion at the 
end of Sect. 5.2. 

Before proceeding to a more detailed calculation of the 
emitted spectrum, we will calculate the slope of the spec- 
trum over a large range of frequencies. It is clear from Eq. 
H15(l how the turn-over frequency drops with the height z 
of the slice of the jet. Since z extends over many orders of 
magnitude from zq up to H, there is a large range of fre- 
quencies v for which part of the jet is optically thin and 
the rest optically thick. The height of the transition 
from the optically thin region to the optically thick one 
depends on the frequency and is given by Eq. H15() if we 
solve for z and call it z^, namely 



y + 4 
K + 3 



(16) 



If, for the moment, we take into account only the op- 
tically thin part of synchrotron emission, to calculate the 
power emitted at frequency v, one has to integrate Eq. 
H13(l from z^ to H and use Eq. H16|l to arrive to the scal- 
ing 



(17) 



The slope ^^^q depends rather weakly on a and for 
3<Q!<5, the slope is ^ 0.2. This result, which will be ver- 
ified by the more detailed calculation of Sect. 5.2, shows 
that the emitted spectra of our model agree very well with 
the flat to inverted radio spectra that are typically ob- 
served from the steady jets of black hole candidates in the 
low/hard state. 

5.2. The radio spectrum 

Knowing the power emitted per unit frequency in each 
slice of the jet (Eq. ^I'S^ ) and the turn-over frequency as 
a function of z (Eq. H15(l ^. the spectrum of the radiation 
emitted by this slice can be approximated by the optically 
thin emission for v > vt that is smoothly connected to the 
optically thick one for v < lyt- Adding the contribution 
of each slice of the jet from zq to H, we have the total 
synchrotron emission in the rest frame of the jet, P{i'). 

For two symmetric jets streaming in opposite direc- 
tions, the total flux that reaches the observer located at 
distance d from the black hole in a direction that makes 
angle with the axis of the approaching jet is (see, e.g., 
Mirabel & Rodriguez 1999) 



47rci2 
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Fig. 3. The emitted synchrotron flux (in mJy) as observed 
at a distance of 1 kpc from our jet is presented for different 
values of the index a and the magnetic field strength at 
the base of the jet Bq. The rest of the parameters are kept 
at their reference values (see Sect. 3). The dotted curve 
corresponds to a = 3 and Bq = 2 ■ 10^ G. The dashed 
curves correspond, in order of increasing flux, to a = 4 
and Bo = 2•10^ 5-105, and 2- 10*^0 respectively. The solid 
curves correspond to a = 5 and Bq = 2-105, 5- 10^, and 2- 
lO^G in order of increasing flux. The spectrum in the radio 
is almost flat with spectral index that depends weakly on a 
(see Eq. (|17|) '). For frequencies above ~ lO^^jJz the whole 
jet is optically thin to synchrotron absorption and has 
slope (1 — a)/2. 

where s is the energy spectral index A log[P(i^)]/A log(i^). 
The Doppler formulae for the approaching, Japr, and the 
receding, ^rcc, jets are given by 



1 



^rcc 



7||(1 — f||/ cosd) 

1 

7||(1 + U||/ cos6') 



(19) 
(20) 



the bulk Lorcntz factor of 



(18) 



Finally, 7|| — 1/y^l — v-^^/c'^ is 
the jet. 

For cos 9 ~ 0.5 and v\\ — 0.8c, the Doppler shifts are 
(5apr — 1, 5rec — 0.4. So, for an intermediate range of view- 
ing angles of our mildly relativistic jet, the contribution 
of the receding jet to the total emitted synchrotron flux is 
rather small (see Eq. I^Hl)). In principle, we should have 
distinguished between the photon frequency in the ob- 
server and in the comoving frame. The two frequencies 
scale with a Doppler factor (5api - Since, Japr — 1 in our 
study, this distinction is not crucial and has been avoided. 

In Fig. Owe plot the synchrotron flux that reaches the 
observer, for a source that is located at distance d — 1 kpc, 
for different values of a and Bq. The rest of the parameters 
are kept at their reference values. Several comments are 
in order in view of this plot: 
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i) The low- frequency part of the spectrum (10^ — 
lO-'^^Hz) has a power-law form with slope given by Eq. 

m- 

ii) The low-frequency cutoff (at ~ 10® Hz) is given by 
the value of the turn-over frequency at the outer edge of 
the jet at height H, while the high frequency cutoff (at 
~ lO^^Hz) is close to the turn-over frequency at the base 
of the jet zq. This break typically lies close to the optical 
region. 

iii) For frequencies >10^^ Hz, the whole volume of the 
jet is optically thin to synchrotron absorption and the 
spectrum has a power law shape with slope — (a — 1)/2 (see 
Eq. 221 )■ The high frequency spectrum steepens further 
as a result of the fact that the energy distribution of the 
electrons extends up to 7max- Equation is accurate up 
to frequency v = {3/2)^^^^h'i^, where i^l = eB / (2nm^c) is 
the Larmor frequency. 

iv) When normalized to a distance of 1 kpc, the ra- 
dio emission of black hole binaries in the low/hard state 
reaches fluxes of up to a few 100 mJy (see Fig. 6 in 
Gallo et al. 2003) . Magnetic field strengths of the order of 
10^ G are required at the base of the jet to explain the 
most luminous (in the radio) sources for 3 < a < 5. 

The part of the spectrum close to the high-frequency 
break of the flat part (close to 10^*^ Hz) is rather approxi- 
mate. Close to the base of the jet the characteristic radius 
of the jet is comparable or even larger than its height, mak- 
ing our analysis inaccurate for this region. Furthermore, 
close to the base of the jet, the electron scattering opti- 
cal depth is high enough that it cannot be neglected. In 
fact, the synchrotron flux emitted close to the base of the 
jet may contribute significantly to the soft photon input 
that is Comptonized by the jet. In the next Section where 
comparison between the spectra calculated from the model 
and the observational data is made, the synchrotron pho- 
tons are also taken into account as a soft photon source 
to be Comptonized in the jet. 

6. Applying the model to XTE J1118+480 and 
Cygnus X-1 

The synchrotron flux, as predicted by our jet model, peaks 
at ~ 10^'^ Hz. Most of the contribution at this frequency 
comes from the lower region of the jet. On the other hand, 
synchrotron photons that arc emitted in this region have 
a significant probability of being scattered by electrons 
once or more times. So, synchrotron emission must also 
be taken into account along with the disk emission as a 
soft photon source. 

For the black-hole binary XTE J1118+480, simultane- 
ous (or nearly simultaneous) observations have been con- 
ducted on multiple occasions at radio, infrared, optical, 
UV, EUV, and X-ray wavelengths (Hynes et al. 2000; 
McClintock et al. 2001b; Frontera et al. 2001), making 
it an ideal source to test the model in the whole elec- 
tromagnetic spectrum. In Fig. 4 the most complete spec- 
tral energy distribution of XTE J1118-I-480 is shown (as- 
sociated with the so-called "epoch 2"), where the radio 



data are from Fender et al. (2001; we do not include the 
observational point at 350 GHz measurement which was 
not done simultaneously with the others) and the infrared 
to X-ray data from McClintock et al. (2001a). Since the 
EUV spectrum depends sensitively on the assumed Nn 
which is not well constrained but probably lies in the range 
1.0-1.3 X lO^^cm-^ (McClintock et al. 2001b), we choose 
to plot the data corrected for A^h = 1.0 x 10^"cm^^ and 
iVn = 1.3 X 1020cm-2. The spectrum of XTE J1118-h480 
has been fitted by an accretion-jet model (Yuan et al. 2005; 
see also Esin et al. 2001), a synchrotron model (Markoff et 
al. 2001) and a thermal-Comptonization model (Frontera 
et al. 2001); for a comparison among different models see 
Chaty et al. (2003). Notice that in the ADAF and corona 
models the presence of a jet is needed to explain the ob- 
served emission in the infrared and radio wavelengths. 

In Fig. 4, the spectrum from our model is also shown. 
The parameters that have been used are Rq = 70rg, a = 
3.5, Tbb = 10 eV and Bq = 10^ G, while the rest of the 
parameters have been kept to their reference values. The 
overall agreement between the calculated and observed 
spectra is very good. Both disk and synchrotron photons 
from the jet contribute an important fraction of the soft 
flux to be Comptonized in the jet. 

It is also clear from Fig. 4 that in this model the opti- 
cally thin synchrotron emission has only a minor contribu- 
tion to the X-ray spectrum in contrast to the jet model of 
Markoff et al. (2001). The main reason for this discrepancy 
is the steeper distribution function of the emitting electron 
that we favor (which results in steeper optically thin syn- 
chrotron spectra) and the higher electron densities in our 
jet model (that result in efficient inverse Compton scat- 
tering in the jet). Recently, Heinz (2004) has shown that, 
adding radiative cooling to the scale invariance formalism 
of Heinz & Sunyaev (2003), which is in accordance with 
the "fundamental plane" of black hole activity (Merloni et 
al. 2003), Compton scattering is favored over synchrotron 
as the emission mechanism of the hard X-rays for the ma- 
jority of these sources. 

Besides its broad-band spectral information, XTE 
J1118+480 has also been extensively observed in the tem- 
poral domain where a large number of interesting proper- 
ties has been revealed for the lightcurves in different en- 
ergy bands. For example, quasi-periodic oscillations have 
been observed with the same frequency in the optical, UV 
and X-rays (Hynes et al. 2003), and the correlation be- 
tween emission at different wavelengths is clear and puz- 
zling (Kanbach et al. 2001; Spruit & Kanbach 2002; Hynes 
et al. 2003). It is interesting (though not within the scope 
of this work) to see whether our jet model is able to re- 
produce these observations. 

We now proceed to apply the model to the December 
1996 RXTE observations of Cygnus X-1 when the source 
was in the low/hard state. For the reference values of 
the parameters, the simulated spectra are in rather good 
agreement with the observed ones (see Fig. O. On the 
other hand, the model underpredicts the emitted flux for 
E<30 keV. 
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Fig. 4. The broad-band spectrum from our model for 
i?o = 70rg, a = 3.5, Ttb = 10 eV, and Bq = 10^ G with 
the rest of the parameters kept at their reference values 
(solid curve) is overplotted with the radio to X-ray spec- 
trum of XTE J1118-I-480. The dashed, dash-dotted and 
dotted curves show the synchrotron, Comptonized and 
disk components respectively. Notice that the EUV spec- 
trum depends sensitively on the assumed Nh and is de 
rived for Nn 

(McClintock et al. 2001b) 



1.0 X 102"cm-2 and iVn = 1.3 x lO^Ocm-^ 



This is the region where the reflection component, con- 
sidered to conic from hard x-rays that are reprocessed or 
reflected by the thin disk, is expected to contribute most. 
In this work, the reflection component has not been in- 
cluded in the calculation. Nevertheless, we have computed 
(using the Monte Carlo simulation for the reference val- 
ues of the parameters) that, assuming that the disk is in- 
finitely thin and that the jet is at an angle of 90 degrees to 
the disk plane, ~ 4% of the X-ray flux hits the disk. If the 
disk is flared or warped (e.g. Dubus et al. 1999) then ~ 8% 
of the X-ray flux reaches the disk (taking an "effective" 
h/R = 0.2, where h is the disk half-thickness and r the 
radial distance) . The fraction of X-rays that hits the disk 
can increase further if the jet is misaligned with respect 
to the outer disk (Maccarone 2002). Thus the reflection 
component can be rather strong and depends sensitively 
on the accretion-jet geometry and the velocity of the jet. 
More detailed modeling is needed to accurately compute 
the strength of this component as predicted by the model. 

7. Conclusion 

The hard X-ray emission in black-hole binaries is usually 
modeled to come from inverse Compton scattering of soft 
photons in a hot corona (e.g. Poutanen & Fabian 1999) 
that lies above and below a thin disk (Shakura & Sunyaev 
1973) or in a two-temperature flow (Shapiro et al. 1976). 
Here we explore an alternative picture where the inverse 
Compton scattering takes place in the jet whose presence, 
whenever the source is in the low/hard state, is well es- 




0.001 O.OI O.I 1 10 
E(keV) 



1000 



Fig. 5. The photon-number spectrum multiplied by for 
the reference values of our jet model (solid curve) is over- 
plotted with the spectrum of Cygnus X-1 as observed by 
RXTE in December 1996. The dashed and dashed-dotted 
curves show the disk-|- Comptonized disk and synchrotron 
self Compton components respectively. 



tablished (Fender 2001; Stirling et al. 2001; Gallo et al. 
2003). 

In Paper I, it was shown that, if the density profile of 
the jet drops inversely proportionaly to distance, power- 
law spectra with photon number index 1.5 - 2, a cutoff at a 
few 100 keV and phase lags almost independent of Fourier 
frequency can be derived for a wide range of the model 
parameters. Furthermore in Paper II, assuming that the 
more energetic electrons lie at small polar distances in 
the jet, both the hardening of the high-frequency power 
spectrum and the narrowing of the autocorrelation func- 
tion with increasing photon energy could be reproduced, 
in agreement with observations of Cygnus X-1. 

Here, we modified this jet model to include a steep 
power-law electron distribution (with 3 < a < 5, that is a 
result of fast inverse Compton -and synchrotron- cooling 
of the electrons) in the jet and checked that all the results 
derived in the X-ray region still hold. We also modeled the 
emission that comes from the jet at longer wavelengths in 
terms of synchrotron emission of the same population of 
electrons and found that the flat to inverted radio spectra, 
observed in black-hole candidates in the low/hard state, 
are a natural outcome of the model. Magnetic fields of 
the order of 10^ — lO^G at the base of the jet are needed 
for the emitted power in the radio to match the typically 
observed one. 

The flat radio spectra have been shown to extend up 
to about the optical wavelengths, with the base of the 
jet emitting mostly in this energy band. This synchrotron 
flux may be strong enough to contribute significantly to 
the soft photon input that is Comptonized in the jet. The 
model has been applied to the broad band spectral energy 
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distribution of J1118-I-480 and Cygnus X-1 showing good 
agreement with observations. 
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